Creation and detection of a mesoscopic gas in a non-local quantum superposition 
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We investigate the scattering of a quantum matter wave soliton on a barrier in a one dimensional 
geometry and we show that it can lead to mesoscopic Schrodinger cat states, where the atomic gas 
is in a coherent superposition of being in the half-space to the left of the barrier and being in the 
half-space to the right of the barrier. We propose an interferometric method to reveal the coherent 
nature of this superposition and we discuss in details the experimental feasibility. 
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It is now possible to control the strength of the atomic 
interaction in a gas, with Feshbach resonances. This has 
allowed the observation of single matter wave bright soli- 
tons with thousands of atoms [l| or a train of solitons 
[3] with 7 Li atoms trapped in a one-dimensional (ID) 
geometry. These solitons are quantum bound states of 
a mesoscopic gas, which opens up fascinating possibili- 
ties: Apart from testing mean field predictions in these 
systems [3J, one can address truly quantum problems, is- 
suing from the quantum nature of the gas center of mass. 

In particular, it was recently proposed to use a Bose- 
Einstein condensate in interferometric experiments to 
test the existence of decoherence mechanisms not pre- 
dicted by usual quantum mechanics and that would show 
up for very massive particles Experiments have suc- 
ceeded in observing interferences with molecules as big as 
fullerenes and there is a need for more massive interfero- 
metric objects [5}. A soliton with a small number of 100 
7 Li atoms has the same mass as Cgo, with appealing new 
features: It does not have internal bound states other 
than its ground state, it can be reversibly dissociated in 
an unbound atomic gas via a Feshbach resonance, and it 
allows the exploration of a new regime, in which the cen- 
ter of mass kinetic energy of the interfering object is of 
the same order as the binding energy of its constituents. 

Furthermore, thanks to the extremely low tempera- 
tures accessible in atomic gases, down to 0.45nK km, and 
the weak decoherence present in these systems [7J, one 
may hope to split the center of mass wavefunction of the 
solitonic gas in two wavepackets that would keep their 
mutual coherence over mesoscopic distances, say a frac- 
tion of a millimeter, much larger than the size of the soli- 
ton. The gas would then have simultaneously non-zero 
probability amplitudes of being in two different spatial 
locations, thus forming a mesoscopic Schrodinger cat in 
reai space. One may then ascertain the presence of a cat 
state by recombining and interfering these two mesoscop- 
ically different quantum states of the gas. This would 
constitute a generalization to many atoms of the one- 
ion experiment of [![. While mesoscopic Schrodinger cat 
states have been reported for radiation fields Q they have 
not been reported yet with ultracold atoms, and atom op- 



tics with a quantum soliton is a promising alternative to 
existing ideas for cat production in these systems [loj ]. 

The dynamics of the center of mass wavepacket during 
the scattering of the soliton on a barrier raises non trivial 
theoretical issues, since the presence of the barrier makes 
the ID many-body problem non integrable via the Bcthc 
Ansatz. We thus construct an approximate effective low- 
energy Hamiltonian for the center of mass of the gas, and 
we derive a rigorous upper bound on the resulting error. 

The starting point is the many-body Hamiltonian in 
ID, for N bosonic particles of mass m interacting via 
the usual contact interaction of coupling constant g, in 
presence of the barrier potential U(x) > 0: 



A' 



H = J2 



2m 



U(xi) 



i<j 



(1) 



This is conveniently rewritten as H — P 2 / (2M)+H- ln +V , 
singling out the kinetic energy of the center of mass 
(M = Nm is the total mass and P the total momen- 
tum of the gas), the so-called internal Hamiltonian H m 
and the sum of the N barrier potentials, V. Without 
a barrier (V = 0) there is full separability between the 
center of mass and the internal variables, so that we split 
the Hilbert space as a tensorial product of center of mass 
and internal variables. H m is diagonalized with the Bethe 
ansatz 11]: For N fixed, its ground state is its single dis- 
crete eigenstate, the quantum soliton \<f>) of energy E (N) 
[13], separated from a continuum of solitonic fragments 
by an energy gap which is minus the chemical potential, 

mg 2 N(N - 1) 



H = E {N - 1) - E (N) 



8h 2 



(2) 



In presence of a barrier, we consider the scattering state 
\ty) of the soliton with an incoming center of mass 
wavevector K > 0. We restrict to a low incoming ki- 
netic energy to have elastic scattering, 

h 2 K 2 

E-Eo^ — KM. (3) 

Far from the barrier, one can then observe only a non- 
fragmented soliton, to the right with the transmission 
amplitude t, to the left with the reflection amplitude r. 
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In this elastic regime, an effective hermitian Hamiltonian 
may be defined, introducing the projector V — IqoM ® 
\<fi)(4>\ acting as the identity on the center of mass and 
projecting the internal state on its ground state, so that 



vm = \<z>) 



(4) 



Far from the barrier, $(X) is simply the center of mass 
wavefunction, X being the center of mass position. The 
so-called VGV formalism, where G is the resolvent of the 
full Hamiltonian [l3j |. then gives the exact equation 

h 2 K 2 



2M 



P 2 

2M +l '< A '» 



sv 



(5) 



The first contribution to the effective potential, in the 
right hand side of ([5]), is the convolution of the barrier 
potential with the internal density profile of the soliton: 



V(X) = (<t>\V\cl>) = 



+ OC 



dx U(X - x)p{x\0) (6) 



where p{x\0) is the mean density of particles in the soli- 
ton knowing that the center of mass is localized in X = 0. 
It was calculated with the Bethe ansatz [14| and is well 
approximated for N 3> 1 by the mean field density profile 
p(x\0) ~ 7V/[4£cosh 2 (a;/2£)], where the mean field soli- 
ton size is £ = H 2 / (m\g\N) . The second contribution in 
([5]) involves virtual transitions to internal excited states: 

Q ;QV\4>) (7) 



SV 



VQ- 



EQ-QHQ' 

where Q = I — V. We shall neglect this contribution but 
not without a justification. From the fact that QHQ > 
Eq + a consequence of the positivity of P 2 /2M and 
V, and of the energy gap of fl"; n , we see in the regime ([3]) 
that the operator — SV is positive and bounded as 



SV < W(X) 



W 2 



V{Xf 



\n\ - h 2 K 2 /2M 



(8) 



When one neglects SV in ([5]), the exact &(X) is replaced 
by <§>q(X), which involves the same incoming wave e lKX , 
but with outgoing waves e %K ' x ' whose transmission and 
reflection amplitudes to and rg are only approximate. We 
have rigorously bounded the resulting errors. We discuss 
here only the experimentally relevant case of an even bar- 
rier U(x) = U(—x). Introducing the "small parameter", 
e = M(<f> \W(X)\$ )/(h 2 K\t \), we have for e < 1/2 the 
following theorem: 



\t-t \ and \r-r \ < 



\to\e 



(9) 



1 - 2e 

It remains to calculate W(X). We have derived from the 
Bethe ansatz the large N asymptotic expression I5j 

W(X) ~ l 2R2 / dx / dyU"(X + x£) 



2M 



— oo J X 



x U"(X + yO 



2 + y-x 



(ev + l)(e~ x + 1)' 



(10) 



In practice, the barrier U(x) is produced with a Gaus- 
sian laser beam, U(x) — Uq exp(— 2x 2 /b 2 ), with a waist 
b much larger than the soliton size £. Then the mean 
potential V(X) is close to NU(X). We shall also assume 
that the incoming kinetic energy h 2 K 2 /2M is about half 
the gap \fi\ ~ ?i 2 /8m^ 2 , so that © is satisfied without 
paying the price of very slow soliton velocities. Then 
Kb ^> 1 and the scattering is in the semi-classical regime 
171 ] . where approximate expressions can be obtained for 
to and rg- A transmission probability 1/2 is predicted to 
be achieved for an incident wavevector Kq such that 



2M 



maxFpO ~ NU . 



(11) 



In the vicinity of K = Kq, the transmission probability 
varies sharply from zero to unity, 

1 



1 + exp[ 



K -K ] 
SK 1 



with SK : 



71V2V 



(12) 



It remains to estimate the bound ([5]). One may take 
U" ~ U"(X) in idDJ), since b > £, so that 



W{X) 



n 2 K 2 

2M 



[u"(x)Y 



^ + 4C(3) 



(13) 



In K = K Q , for e <C 1, a semi-classical calculation gives 



\t-t \ 



< 



mm 3 

N 1 / 2 



In (Nb 2 /£ 2 



(14) 



a quantity checked to be -C 1 in what follows. 

We now study the experimental feasibility. An axial 
Gaussian laser beam confines N ~ 100 atoms of 7 Li in 
the y — z plane, with a resulting transverse harmonic 
oscillator length = (h/muj^) 1 ^ 2 ~ 0.54/xm, where 
wj_ — 27r x 4.8KHz is the transverse oscillation frequency. 
In this optical wave guide, the interacting gas has a one 
dimensional character if 2£ ^> a±. In order to make 
cooling of the gas not too challenging, we take a not too 
large soliton length £ ~ 0.9/im; the resulting 3D scat- 
tering length, a ~ —a 2 L /(2N^) ~ — 1.72nm is in the in- 
terval of values (—00, — 1.5nm) accessible with the Fesh- 
bach resonance [l[. Initially the gas is also harmonically 
trapped along x with an oscillation frequency lo. The gas 
is assumed to be cooled to the temperature T = 0.45nK 
Q. This axial trap is so weak that it very weakly af- 
fects the internal solitonic variables, fiw < |^|/10, but it 
is strong enough that the center of mass of the gas, still 
separable in a harmonic trap, has a negligible probability 
cxp(— fkij/ksT) < 1/10 to be in an excited state. These 
two constraints impose the weak value u> ~ 2ir x 23.5Hz. 
They also imply |/x|/fesT ~ 25, so that the internal vari- 
ables of the soliton are frozen in their ground state. 

At t — 0, the gas is launched with a total momentum 
hKo such that 



2M 



2 



16m£ 2 



(15) 
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The corresponding velocity is HKq/M ~ 0.37mm/s. Si- 
multaneously the axial trap is switched off, to free the 
center of mass of the gas, with an initial wavepackct 



§(X) OC e ^oX e -(X-X f(^K) 



(16) 



With a sudden opening of the axial trap [18j |. 
h 2 (AK) 2 /2M = hLu/4 and AK/K Q ~ 0.22 < 1: This 
wavepacket is quasi-monochromatic. Smaller values of 
AK may be obtained by a more clever opening procedure 
of the trap, within times ~ I /to 19]. The wavepacket is 
then scattered on a broad Gaussian barrier centered in 
x = (here X < 0), a beam-splitter, created by a laser 
beam of waist b = 5£ 3> £ and of intensity adjusted to 
satisfy the half-transmission probability condition (jlljl . 
In any realistic case, AK remains much larger than SK, 
so the wavepacket experiences a mere filtering in Fourier 
space, the components with K > Kq being transmitted 
and the ones with K < K being reflected :20| . As a 
consequence, the wavepacket also splits in real space in 
a transmitted part and a reflected part, that nicely sep- 
arate since their mean velocity exceeds their spreading 
velocity: A mesoscopic Schrodinger cat is born. 

How to prove this experimentally ? The first step is to 
check the absence of fragmentation: A photo of the gas 
by absorption imaging should show, for any realization 
of the experiment, that all the particles are clustered in 
a single lump of size £, randomly situated to the left 
or to the right of the beam-splitter. The second step 
is to check that the two wavepackets are coherent, by 
recombining them and looking for interference fringes, 
with a fringe spacing tt/Kq. The recombination of the 
two wavepackets is obtained by their total reflection on 
mirrors, produced by two Gaussian laser beams centered 
in x — L/2 and x = —L/2, L ^> 1/AK, with the same 
waist as the beam splitter but with a higher intensity 
(say, twice as high). The reflected wavepackets interfere 
around x — 0, the beam splitter being switched off 21 1. 

We have studied the proposed experiment by a numer- 
ical solution of Schrodinger's equation for the center of 
mass wavefunction, with the initial condition (fT6|) and 
the same approximate effective Hamiltonian as in scat- 
tering theory: 



ihdt$(X,t) = 



-M 



2M 



x 



V{X,t) 



9(X,t) (17) 



The center of mass probability distribution |$(X, t)\ 2 is 
plotted at key times in FigfT] To quantify the contrast of 
the interference fringes, we also plotted the modulus of its 
Fourier transform, s{Q,t) = f+™ dX e' l( ^ x \^{X, t)\ 2 . 
When the two wavepackets overlap, sharp peaks in \s(Q)\ 
indeed form in Q ~ ±2Kq, with a contrast |s| ~ 0.32. 
This is a high value, as the ideal case of two overlapping 



plane waves 3>(X) oc e 



KnX 



iK X 



;ives 1/2. 



The high contrast interference fringes in FigfT] are how- 
ever for the center of mass probability distribution, not 
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FIG. 1: Evolution of the center of mass wavefunction of a 
solitonic gas with = 99 atoms, by integration of (|17|l for 
the initial condition JTSJ, with AK ~ 0.093^0, X =_-15f. 
Left panel: |"l>(X,t)| 2 (solid lines), effective potential V(X,t) 
(dashed lines). Right panel: modulus of the Fourier transform 
s(Q,t) of \${X, t)\ 2 ; the vertical lines are in Q = ±2K . The 
time is in units of ML/hKo- Att — 0.63 the state is a cat. At 
t — 1.31 the two reflected cat components strongly interfere, 
two narrow peaks of height 0.315 emerge on |s(Q)| in Q = 
±2Ko- Maximal interference occurs at t = 1.36. The last line 
is an average over a Poisson distribution for N, with a mean 
value N — 99; the peak height in |s(Q)| is reduced to 0.062. 



for the atomic density, which raises the question of their 
observability by usual fluorescence imaging. The mean 
atomic density p(x) is the convolution of |$(A)| 2 with 
the internal soliton density p(x\0); since the soliton size 
£ is as large as the fringe spacing tt/Kq, one finds that the 
contrast of the fringes in p{x) is several orders of mag- 
nitude smaller than in |<I>(A)| 2 . This problem can be 
solved by increasing, just before imaging, the intensity of 
the laser producing the transverse trapping by a factor 
about 21, which reduces the transverse harmonic oscilla- 
tor length to d± = 0.25^m and brings the soliton close 
to its collapse threshold N\a\/d± ^ 0.67 Further- 
more fluorescence imaging can be optimized to measure 
directly the quantity |s(2i£"o)|, by exciting the gas with a 
laser standing wave along x, produced by the superposi- 
tion of two laser waves of wavevectors k± = (±k x , k y , 0) 
at some angle with the x axis such that k+ — k- = 2Kqc x . 
The resulting fluorescence rate in direction n per unit 
of solid angle is given in the Born approximation by 



dT/dfl oc 



Ei=i e 



with validity conditions 



discussed in [23j. Here e(r ) is the laser electric field. The 
emission rate Tq of photons in the solid angle of the 
detection lens is an oscillating function of the location of 
the antinodes of the laser standing wave with respect to 
the interference pattern in |<f>(A)| 2 , with a contrast 



rmax T^mm 

1 n ~ 1 n 

rmax i "pmin 

1 n "+" 1 a 



= \s{2K )\S in (Cl). 



(18) 
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The reduction factor 5i n (fi) is a function of the 3D static 
structure factor of the soliton for fixed center of mass 
position, that we approximate with the 3D mean field 
theory. By using a lens of optical axis along k + + fc_ 
with a numerical aperture 0.4, one finds the remarkably 
hig h value S- m (Q) = 0.84, thanks to a superradiant effect 
23j, which also concentrates 16% of the fluorescence in 
the 4% solid angle fraction collected by the lens. 

It remains to check that decoherence is negligible dur- 
ing the transit time t% Iaxls — ML/TlKq < 200ms of the 
cat state in the interferometer. In cold atom experi- 
ments, the main source of decoherence is particle losses: 
A single loss event would destroy the cat, since it "mea- 
sures" the positions of one or several atoms and local- 
izes the center of mass of the gas within the soliton 
size £. The usual loss rate formula for m-body loss is 
dN/dt = —K m J d 3 rn m (f); here one should take for n 
the 3D density profile for a fixed center of mass position, 
that we approximate with the mean field theory. For one- 
body losses due to collisions with the background gas, one 
should have a loss probability KxNttmns < 1/10, which 
imposes the reasonable lifetime K^ 1 > 200s. For three- 
body losses due to formation of deeply bound dimers, 
the loss constant K$ for 7 Li at the considered magnetic 
field B is not known. Since |a| is smaller than the Van 
der Waals length 3nm, as it is for B — 0, we use the 
B = prediction of [24|, applying the factor 6 reduction 
for a condensate, K3 s» 3 x 10 _41 m 6 /s, which leads to a 
negligible loss event probability ||dJV / dt\t tl ans ~ 0.03. 

In present experiments the number of atoms N fluctu- 
ates from one realization to the other around the desired 
mean value N. Since the launch velocity HKq/M is fixed, 
Kq is proportional to N and also fluctuates [25|. A first 
side effect is that the half-transmission probability con- 
dition may be violated away from TV = TV; fortunately 
this is not the case for a broad barrier b ^> £, since both 
terms of (fTTj) are proportional to N. A second side effect 
is that the fringe spacing tt/Kq will fluctuate, which will 
blur the fringes. A simple way to estimate this is to as- 
sume that |$(A)| 2 cx \e lK ° x +e -K x ] 2 e -x 2 /2^ at the 

overlap time. Averaging over a Poisson distribution for 
N with cr and K Q /N fixed leads to, for \X\ < nN/K : 



e -X 2 /2a 2 

(2 7 r) 1 /V 



,-X 2 /2cr 2 



cos(2K X) 



The fringes persist around the origin over a distance cr c = 
N 1 / 2 /(2K a ) = V2£. \s{2K )\ is then reduced by a fac- 
tor ct c /(ct 2 + cr 2 ) 1 / 2 . Estimating a from FigfJJleads to a 
reduction factor 5 close to the numerical one (see FigQ}. 

In conclusion, we propose to produce a coherently bilo- 
calized gas by scattering an atomic quantum soliton on a 
barrier. We have performed a detailed analysis of this 
idea, which raises challenging experimental aspects of 
preparation and detection, but also non trivial theoreti- 
cal aspects since this is a many-body problem. We find 



that a gas with N ~ 100 7 Li atoms can be prepared in a 
coherent superposition of being at two different locations 
separated by ~ 100/im, and that this can be proved by 
an interferometric measurement. 
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